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Abstract 

The Marginal Fermi Liquid (MFL) hypothesis proposed by Varma 
et.al is derived from a consistent treatment of local spin fluctuations 
in the Falicov-Kimball (FK) Model. Within an infinite-dimensional 
mean- field approach, which is exact for this model in d = 00, a 
marginal Fermi liquid susceptibility and single particle self-energy are 
obtained near half-filling, and a Fermi liquid phase is recovered suf- 
ficiently away from half-filling, in agreement with indications from 
experiments which probe the normal state of cuprates as a function 
of doping. 
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The discovery of high-T c superconductors in Cu — O based compounds 
has led to an upsurge of theoretical work concerning the unusual normal state 
properties of these materials, which appear not to conform to the framework 
of the Landau Fermi liquid theory [1, 2]. A way of unifying the diverse 
anomalies observed in experiment was proposed by Varma and coworkers [2], 
who suggested a phenomenological ansatz for the spectrum of charge and 
spin fluctuations. For low frequencies w C ^ this marginal- Fermi-liquid 
(MFL) ansatz is 



where oo c is a cut-off energy. The s.p. self-energy cr(u) ~ ulnuo ± i\ou\, as a 
consequence of (1), and this reconciles the unusual normal state anomalies 
with the existence of the Luttinger Fermi surface. 

As emphasized in [2], the singularities in (1) are in the frequency depen- 
dence; the momentum dependence is assumed smooth. In principle, an exact 
solution of certain strongly correlated models in d — oo, where local fluctu- 
ations are treated exactly, should lead to the above spectrum. Furthermore, 
if these singularities do not depend on any special symmetries which are lost 
in the lattice problem, they are likely to survive in the lattice problem. 

Varma et. al [2, 3] have solved the multiband Hubbard model within the 
impurity approximation to obtain the MFL form for the local susceptibilities. 
However, the MFL theory has not been able to account for the T and doping 
(x) dependence of the Hall constant R H . Recently, Mahesh et. al [4] have 
computed the T and x dependence of Rh by numerical diagonalization of the 
one-band Hubbard model on finite-sized clusters. They were able to account 
for the anomalous behavior of Rh- However, a proper description of all the 
anomalous features has not been possible, and the extension of the Luttinger 
liquid [5] concept of Anderson to two dimensions is not clear, unless the small 
momentum forward scattering couplings become singular [6]. 

In this letter, we show that a consistent treatment of local fluctuations 
in the doped, large U Falicov-Kimball model leads to the Eqn. (1). We 
utilize the exact solution of the infinite dimensional effective Falicov-Kimball 
(FKM) model [7]. An explicit analytical calculation of the dynamical spin 
susceptibility Xo-(^) is consistent with the MFL form. 
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In this paper, the Falicov-Kimball model in 2d in the large U limit, 

H = -t ( c l c j + h - c ) + U Y1 n ic n id ~ ^^{n ic + n id ) (2) 

<ij> i i 

is proposed as an effective model capable of describing the anomalous prop- 
erties of oxide superconductors in their normal state. In this model, t and 
U should be understood as effective parameters which are determined by 
comparison of the low-energy spectra of (2) with that of a more realistic 
three-band model [2,8]. The proposed FKM bears some similarity to the 
effective model (eqs (12)-(14) of ref. [2]). However, since the authors of ref. 
[2] solve an impurity model, they require fine-tuning of parameters to reach 
the critical point. Since we perform a lattice calculation exact in d = oo, the 
critical behavior survives for a finite range of filling, as the authors of ref. 
[2] anticipate. In this FKM [7], the d holes are immobile. This means that 
[n id ,H] = OVi; hence, the model has an exact local U(l) symmetry. We no- 
tice that the model eqn (2) is different from the usual Hubbard model, which 
has a global U(l) symmetry associated with total fermion number conserva- 
tion. As we shall see, it is this local symmetry which leads to the breakdown 
of Fermi liquid theory in our approach. 

Since we are interested in the nontrivial local dynamics, we consider an 
auxiliary impurity model in which the d hole does not hybridize with the 
"conduction electrons". This impurity model hamiltonian is the d — oo 
counterpart of Eqn. (3), and reads 

H = J2 £ kc\ck + t'$2(e %k ' Ri c\c k + h.c) + Un ic n id -^n(n ic + n id ) (3) 

k k i 

where i represents the impurity site, and k labels the "conduction electrons". 
We are interested in the non-trivial local dynamics; hence, we compute the 
local s.p and two-particle propagators exactly. To compute the s.p Green 
function, we start with an equation of motion for it. 

(iu e + njG^iuje) = + *fc G w( iw <) + u (( n idCi] c\)) iu!e 

271 k 

(iuj e + /J, - U)((n id c i ]c ] i ))i Wt = < U ' d > +y2t k {{n id c k ; cf))^ 

2n k 
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(iu) t - e k )G c ki (iu e ) = t k G c u {iuJt) 
(iu e - e k ){{n id c k ;cl)) iu!e = t k ((n id a; c}))i 



(4) 



solving for G^(iu^) yields (for a Lorentzian unperturbed DOS with a half- 
width A, Po (z) = (A/n)(z 2 + A 2 )- 2 ) 
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2tt 

with the self-energy 



1 -n d 



n id 



iui + ji + iAsgriUi iui + /i — £7 + iAsgnun 

U 2 n d {l - n d ) 



T, c {iu t ) = -ji + Un d + - 



iuj£ + ji — C/(l — n d ) + iAsgncvg 



Also, it is easily seen that 



1 



27r iu^ + // — [/ + iAsgnue 
The s.p and the two-particle local spectral densities are 



A 



1 -n d 



+ 



and 



2vr 2 L(^ + /i) 2 + A 2 (^ + - t/) 2 + A 2 _ 
A 



P {2) (^) 



(5a) 



(5b) 



(6) 



(7) 



(8) 



2tt 2 (jut + fi-U) 2 + A 2 

^From Eqns. (7) and (8), it is clear that the low-energy spectrum is a super- 
position of s.p and two-particle states. It is precisely the resonant scattering 
between these states that leads to the breakdown of FLT. 

To proceed with the derivation of the MFL spectrum, we compute the 
NMR relaxation rate, which is related to the low-frequency dynamical local 
spin susceptibility via 

1 = -TlimE^,for W «T 



= -Ex"(?»- forT«c. 



(9) 
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where "x!'{oj) is the imaginary part of the dynamical spin susceptibility. We 
also have [91 
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A 2 



2hN 
A 2 



E r°°{T T [St{r)Sr m e^dr 

• J — oo 



uj->0 



2nHAN z - 



(10) 



so that the task reduces to calculating (exactly within d = oo) the lo- 
cal, transverse spin correlation function. This can be computed easily once 
p( 2 \iut) is known, by using the identity 



(S?Sr) 



n d 



(ri ic n id ) 



(11) 



But D = (n ic nid), the average number of doubly occupied sites is com- 
puted directly from the two-particle spectral density via a Matsubara sums 

1 



(nicUid) 



P 



Ep (2) (^)e- W(H) 



Direct evaluation by use of Eqn. (9) yields 



{ni C n id ) = — 



7T L2 



- + tan" 1 ' p 



7T 



+ tanh|/5 



A 
H-U 



Hence from Eqns. (12) and (14), we get 



n d , _i (Un d 
— tan — - — 
7i V A 

JJn, 



k B T < A 
, k B T^>A 

k B T<^A 



= ^tanh(^) , A; B T»A 



(12) 



(13) 



(14) 



where we used the exact relation /j, — U(l — nd). Comparison of Eqn. (10) 
with Eqn. (11) (after substitution of (15)) leads to 



X» = -^(O^tan- 1 (^j (^) , u«T 



5 



A 2 

= -—p{Q)n d , PU > 1 , T < uj 



(15) 



This is precisely of the MFL form. A second-order perturbative calcu- 
lation gives S(cj) ~ uoinuo — i\uo\. Thus, a proper treatment of local spin 
fluctuations leads to the MFL spectrum for the spin susceptibility as well 
as the self-energy. It is important to notice that the MFL ansatz is a state- 
ment about both the s.p self-energy as well as the susceptibilities, and so 
approaches which compute only the self-energy of the correct form [10] are 
inadequate. 

Earlier attempts, besides those of [2, 3] have attempted to derive the MFL 
spectrum by invoking a negative U HM [11], or the Holstein model [12]. The 
relevance of these models is questionable, since there are no antiferromagnetic 
insulating or spin fluctuation dominated strange metallic phases in these 
models. Our approach, which starts from the 2d, Falicov-Kimball model 
with large U, explicitly takes AFM local spin fluctuations into account. In 
fact, noticing that the FKM in d = oo is the recoilless x-ray edge problem 
[13], the local "excitonic", or the transverse spin-spin correlation function is 
divergent at low energy near n — 1, 

POO 

**» = / dre^{T T [S + {r)S-m) 
Jo 

~ \<v\- p (16) 

with P — ^ — (J^j and 5 = tan -1 (^Q, the s-wave phase shift at /i. This 
leads to a soft, local spin fluctuation mode at low energy, and it is precisely 
the coupling of the s.p part of the spectrum to these soft modes which leads 
to the MFL ansatz. 

The MFL ansatz reconciles the normal state anomalies observed in cuprates 
with the existence of the Luttinger Fermi surface. We have shown that the 
MFL ansatz can be derived from a Falicov-Kimball model. Large U is cru- 
cial to the derivation, as is the filling factor. The MFL state is unstable to a 
Fermi liquid (FL) phase for n < n c = 1 - ^cor 1 (^) [13]. With U = 2 A, for 
e.g., n c = 0.75, i.e. the doping concentration of holes x = 0.25, close to the 
experimental x c value beyond which the anomalous behavior is suppressed 
[14]. At lower densities, effects of disorder will introduce qualitatively new 
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features. A non-perturbative treatment including effects of disorder will be 
reported elsewhere. 

We have not dealt with the superconducting phase in this paper; however, 
a few remarks can be made. Due to large U, the on site pairing amplitude 
is severely suppressed, (c\d\) = (actually, because the Hamiltonian has 
an exact local U(l) symmetry, this is rigorously true, by Elitzur's theorem). 
This implies, upon Fourier transforming, that 

£<ct(*)<*t(-£)> = o (17) 

k 

As pointed out by various authors [15, 16], the pair wave-function should 
have lines or points on the FS at which it is zero, and hence the pairing 
is not of the BCS variety. The gap function has nodes at the same points 
where the pair wave function does. It is not possible, from the above, to 
specify whether the symmetry of the SC state is of the extended s-wave or 
the d-wave type. This is beyond the scope of the present work. 

In conclusion, we have shown that an exact treatment of the dynamical 
spin fluctuations in the large U, one band Falicov-Kimball model leads to the 
MFL phase near half-fillings and a FL phase sufficiently away from n — 1. 
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